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Abstract 



One loop corrections to the domain-wall quark propagator are calculated 
in massless QCD. It is shown that no additative counter term to the current 
quark mass is generated in this theory, and the wave function renormalization 
factor of the massless quark is explicitly evaluated. We also show that an 
analysis with a simple mean-field approximation can explain properties of the 
massless quark in numerical simulations of QCD with domain-wall quarks. 
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I. INTRODUCTION 



The formulation of the lattice fermion in QCD with the chiral symmetry is one of the 
most fascinating problems theoretically and practically. Although both Wilson and Kogut- 
Susskind (KS) fermion formulations have been popularly used for the lattice QCD simula- 
tions, some disadvantages remain in these formulations: In the Wilson fermion formulation 
the quark mass has additive quantum correction and the chiral limit is reached only by the 
fine tuning of the mass parameter. As a general rule we have to take continuum limit tuning 
the mass appropriately in order to simulate massless QCD. In the KS fermion formulation 
the number of flavors is restricted and the original flavor symmetry is broken explicitly to 
some residual one. 

The domain-wall fermion formulation, which was originally proposed to define lattice 
chiral gauge theories [0], has been applied to the lattice QCD 0. This formulation is 
expected to have great advantage over the previous two formulations: An advantage over 
the KS fermion is that the number of flavor is not fixed. This is manifest from its definition. 
The other advantage over the Wilson fermion is that mass renormalization is multiplicative 
(m e s = Z m mtree)- In other words, if a massless mode exists at the tree level it is stable 
against the quantum correction. This property is not a trivial one, but only an intuitive 
discussion on it has been given so far On the other hand, the recent numerical simulation 
suggests that the stability of the zero mode holds even non-perturbatively [|J. Therefore an 
analytical understanding of the domain-wall QCD is now needed. The aim of this paper is 
to confirm the stability of the massless mode by the lattice perturbation theory and to give 
explicitly the wave function renormalization of the quark field. 

This paper is organized as follows. In section 2 we will give basic tools for the perturbative 
calculation with the domain-wall fermion. It is enough to present only fermion propagator 
because other Feynman rules of gauge interaction and gauge propagator are exactly identical 
to that of the ordinary Wilson fermion. In section 3 we calculate one loop corrections to 
the fermion propagator. Section 4 is the main part of this paper, where we discuss the 
renormalization of the zero mode or massless quark field. We take the diagonal basis of 
the mass matrix of the domain-wall fermion and see that the zero mode is stable against 
the one loop correction. The wave function renormalization factor of the massless quark 
field is also given explicitly. Section 5 is devoted to the mean field analysis. We show that 
properties of the zero mode observed in the numerical simulation || are well explained in 
this approximation. In section 6 we give our conclusion and discussion. In appendices some 
derivations of formulae used in the text are presented. 

In this paper we set the lattice spacing a = 1 and take the SU(N C ) gauge group with 

N 2 - 1 

the gauge coupling constant g and the second Casimir C 2 = — ^— - — . We set iV c = 3 in the 
numerical calculations. 

II. PERTURBATION THEORY WITH DOMAIN- WALL FERMION 
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A. Action 



We adopt the domain-wall fermion of Shamir type to describe massless quarks. The 
domain-wall fermion is a variant of the Wilson fermion with sufficiently many flavors and 
special form of the mass matrix. Although it is also interpreted as a five dimensional Wilson 
fermion QTJ , we prefer to treat it as the multi-flavor system 

In this point of view only difference from the Wilson fermion action is the fermion bilinear 
term. If we separate the QCD action for lattice perturbation theory into fermion and gauge 
parts, 

S 'S'fermion "Sgauge Sqp -\- Spp -\- ^measure; (1) 

the lattice gauge action S gauge , the gauge fixing and the FP-ghost term Sqf + Spp, the in- 
variant measure term S meaBwe and the gauge-fermion interaction terms in Sfermion are exactly 
same as those in the ordinary Wilson fermion perturbation theory |||6| with many flavors. 
The domain-wall fermion action Sfermion is written as 

Sfermion = ^ J> m , a ilti D n) m> n ^n,s + ^s^m'n P +^n,t + ^m.^m.n^-V'n,* 
n,m 

+ m q ij) m s (8 mjn 8 StNs 8 t ,iP+ + 8 m!n 8 Stl 5 t ,N s P~) (2) 

where m,n is four dimensional space index, s, t = 1,---,N S is the flavor index. Here the 
Dirac operator is given by 



, 2 



and mass matrix W s t is defined as 

W n,m;s,t = £*fcl,A,m ~ W n . m 5 s j (4) 

where 

y / \ 

Wn,m = (1 — M)8 n:m + — [Un,ii8n+(i,m + ^m,/i^"--A: m _ ^^n,m) (5) 

Z I* 

is a sum of the Dirac mass term and the Wilson term, which contain gauge fields at this 
stage, and r is the Wilson parameter, which we set r = — 1. The parameter m q is the current 
quarks mass, but in this paper we only treat the massless QCD taking m q = 0. P± is a 
projection operator defined by 

P± = ^ (6) 

In our domain-wall fermion action (Q) we have Dirac mass M besides the current quark 
mass m q . Here we have to notice that M is not the physical quark mass but it is rather 
an unphysical mass of the cutoff order (1/a) like Wilson term. As will be mentioned later 
M has an important role as a parameter of the theory: choosing a suitable value for M we 
have a massless fermion mode for the vanishing current quark mass (m q = 0). 
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In order to see the massless fermion mode it is more convenient to be in the momentum 
representation and pull out the bilinear term. The fermion action in the momentum space 
is written as 



'fermion 



d A p 



£^sin PA1 + W + (p) s , t P + + W-(p) a , t P- 



int. 



where the mass matrix has the following form, 



W + {p) 



s.t 



S a+1 , t - W(p)5 s j 
f-W{p) 1 

-W{p) 



-W(p) J 



V 

W-(p) a , t = Ss-^t - W(p)S a , t 
f-W(p) 

1 -W(p) 



V l -W(p)J 

W(p) — 1 — M — r 2(1 - cosp M ). 



(7) 



(8) 



(9) 



(10) 



The gauge interaction term Si nt . is identical to that of the Wilson fermion perturbation 
theory with N s flavors. 

As will be discussed in appendix C, in spite of the presence of the Dirac mass M this 
fermion system has one massless fermion mode and N s — 1 excited modes with the mass of 
cut-off order by virtue of this mass matrix form, provided that \W(p ~ 0)| < 1 is satisfied by 
a suitable choice of the Dirac mass M. Here we take the momentum region p M ~ to see the 
zero mode with physical momenta. At the momentum p^ ~ n, where the doubler emerges 
in the naive fermion formulation, the parameter condition is not satisfied < — • vr) ] > 1), 

so that all N s fermion modes have mass of the cut-off order. 



B. Fermion Propagator 

In the next section we will calculate the one loop correction to the fermion propagator. In 
this subsection we set up the lattice Feynman rules for domain-wall fermion with vanishing 
current quark mass (m q = 0). 

As is discussed in the previous subsection, the domain-wall fermion action is almost same 
as that of the ordinary Wilson fermion's one with N s flavors. The peculiar feature of the 
domain-wall fermion is the form of the fermion propagator, which is given by 

s F (p) 9 ,t = kjv + w + ( P )p + + w-( P )p] ] (ii) 



s,t 



where p^ = sin(p^). The explicit form is written as 
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S F (p)s,t = [H7„IV + W-)G R (s,t)P + + (-17^ + ^+)G L ( S ,t)P_] st (12) 

where 

1 \ 



G R (s,t) = - 



St 



p 2 + w+w- 

G°{s -t) + A ++ e a(s+t) + A + _e ai - s - t] + A_ + e a{ ~ s+t) + A__e a{ - S ' l) (13) 



G o (s — t) = A ( e °C*H*-t|) + e -«(*.-l«-*l)) (14) 
A ++ \ _ A ((l-We- a )(e- 2aN °-l)\ g) 

(16) 



/ e a7V - (1 - WV*) - e~ a7V - (1 - We~ a ) V W(e° - e 



— a \ 



A__ J e aN °(l - We a ) - e- aN °(l - We~ a ) \ We a )(l - e 2aN °) 



and 

G L (s,t) 



St 



\p 2 + w-w+ 

G°(s -t) + B ++ e a(s+t) + B + _e a{s - t] + B^ + e a{ - s+t) + B^e 01 ^ 8 ^ (17) 



B ++ \ A fe- a {e- a - W){e~ 2aN ° - 1 



B_+ J e aN °(l - We a ) - e- aN °(l - We~ a ) V W(e a - e 

B + \ A ( w(e a -e- Q 



e aN °(l- We a ) -e- aN *(l- We~ a ) V e a (e a - W)(l - e 2aN °) 
Here a and A are defined as 



(18) 
(19) 



cosher I±i^ (20) 
2W(p) V ; 

sinha = ^^(1 - W 2 ) 2 + 2(1 + W 2 ) $>n 2 /v + QT sin 2 ^) 2 (21) 
1 l 

^4 = C22) 

sinh a 2 sinh(aiV s ) V ; 

Note that the argument p of W and a is suppressed throughout this paper unless necessary. 
Since this fermion propagator is invariant under a — > —a, we take the a > without loss of 
generality. G R and are also symmetric in (s,t). See appendix B for the derivation. In 
the one- loop calculation we use the above propagator in the N s — > oo limit. 



III. ONE LOOP CALCULATION 
A. Diagrams 

In this section we calculate the one loop correction to the fermion propagator, which is 
;iven by two contributions £ tadpole (p) _|_ ^haif-arcie^*^ from diagrams in Fig. |l[ 
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The 1PI fermion 2-point vertex function is given by 



Vi-ioop(p) s,t 



i-ff, sinp M + W + (p)P + + W (p)P_ - S(p) 



s.t 



(23) 



with 



£(p) = S tadpole (p) + S 



half —circle 



(P)- 



(24) 



In order to investigate the massless mode of Ii_ loop (p) Sjt in the p^ — > limit, we need 
only the first few terms in the p^ expansion. Since the only dimensionful quantity is the 
external momentum p^ in our calculation, the higher order terms in the p^ expansion are 
also higher order in a. 



B. Contribution from Tadpole Diagram 

The contribution from the tadpole diagram is written as 

S tadp ° lc = \g 2 C, £(*7, sinp M - r cosp M ) £ T^T^lA* (25) 



(27r) 4 4sin 2 | 

! rC 2 T(^ + 2),l. l + 0(o) (2G) 



.-v ' 1 



where T is the tadpole loop integral 

dH 1 



-7T (2vr) 4 4sin 



2 { = 0.154612. (27) 



The first term in Eq. (|26"D is finite, and the second term linearly diverges in the limit a — > 0. 
We see that £ tadpole is diagonal in flavor space, and its effect is to modify the mass parameter 
.\/ — M M - 2g 2 C 2 T. 



C. Contribution from Half Circle Diagram 

The contribution from the half circle diagram in the Feynman gauge 

/•7T rj^l ( 1 1 ~| 

ss if - circie = £ ^ ^ { ~ wTa) p cos 2 {l " + ^ ~ ir sin 2 (/m + 

x 5 F (0^X (-^T a )|7 M cos^(/ / ,+^)-zrsin^(/ M +p M )| x ^-j^ 2 ( 28 ) 

cannot be calculated analytically because of its complicated dependence on the flavor indices 
s,t in the fermion propagator. Here p^ is defined as p^ = 2sin(p M /2). 

It is easily seen that the loop integral of (p8[) has infra-red divergence. As is in the 
ordinary lattice perturbation theory the infra-red divergence can be written in an analytic 
form. To do this we separate E^f f_circle as follows: 



6 



-ihalf — circle / 



(p) = £^-(p) + £«f-(p) 



(29) 



where 

E]5-(p) = E5f- drtto (p)-E^ t -(p) 
and S™ nt -(p) is introduced to extract the infra-red divergence: 



(30) 



^\p)=2g 2 C 2 J 



dH -if{C + P+ + C_P_ 



(27T) 



6(n 2 -l 2 ) 



with (C+) S)t = (1 -w 2 )^ + '~ 2 , (C_) Sji = (l-^)^" 5 -', and w = W(0). In order to have 
zero modes with the physical momentum, wq should be in the region Wq < 1. This leads to 
the condition of M that < M < 2. Since S^ a j: (p) is infra-red finite in p —>■ limit, we can 
evaluate it in the p expansion: 



nylat. 



(31) 



The logarithmically divergent part S cont (p) can be calculated analytically, while a linearly 
divergent and finite terms (the first and the second terms in eq.(^l|)) have to be evaluated 
by numerical integrations of loop momenta. After a little algebra we have 



ihalf — circle 
J s,t 



ip 1 II+P+ + I- t P\ + M+ P + + M~ t P_ 



where and M ± are given by 



1% ~ Aog( S ' t) + ^finite ( S ' 



16tt 2 - 



9 2 C 2 (C ± ) s , t ln(7r 2 ) + --lnp ; 



(32) 

(33) 
(34) 



4nteM) = 9 2 C 2 J lY.( c ™UW-G R + W + G L )(s,t) + sm 2 l,(G L + G R )(s,t) 



(2vr) 4 /2 



+ $>n 2 Z„/2G fl (s,t) 



(W-G fl + W + G L ) (s, t) + 2(^ cos 2 Z„/2 - 2 cos 2 lj2)G L {s, t) 

rf 4 / 1 



-^ 2 c 2 (c + ) S)t I 



(2vr) 4 (Z 2 ^ 



(tt 2 - f ) 



(35) 



WM) = s 2 g 2 J 



dH 1 



(2vr) 4 p 



£ (cos /^(W-Gb + W + G L ) (s, t) + sin 2 Z M (G L + G R ) (s, t) 



y, sin 2 1^ 



2\2 



Kt = 9 2 C2 J 



/. 4(Z 

^sin 2 /,/2G L ( S ,t) 

V 

dH 1 



(W-G* + W + G L ) (s, t) + 2(£ cos 2 l v /2 - 2 cos 2 lj2)G R {s, t) 

V 

dH 1 



-g 2 C 2 {CJ) s>t J 



(vr 2 - Z 2 ) 



(2tt) 4 Z 2 ^ 



(2tt) 4 (Z 2 ) 2 

cos 2 Z M /2(^+G L ) (s, t) - sin 2 Z / ,/2(^"G fi ) {s, t) 



(36) 
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+ ^sin 2 l^G L + G R )(s,t) 



(37) 



M, 



2 f d A l 1 _ 



cos 2 / M /2(^-G R ) (s, t) - sin 2 IJ2{W + G L ) (s, t) 



+ - sin 2 I^Gl + Gr)^,^ 



(38) 



By the dimensional counting 1^ has In a 2 divergence and constant terms in a, and M ± has 
1/a linear divergence when lattice spacing a is introduced explicitly. Although M ± may 
have In a 2 divergence naively, it is canceled by the algebraic relation 



(W + (p — 0)) Sjt (l — Mf — 
(W-(p = 0)) 8>t (l-M)- t = 0. 



(39) 
(40) 



The logarithmic divergence In a 2 in I^ 1 is given analytically. As we can see from the form 
of (C + ) s>t = (1 — w%)wQ +t ~ 2 , l£ is localized in the boundary (s, t) = (1, 1). This is because 
the logarithmic divergence comes from the effect of massless fermion mode which is localized 
in the boundary. The other one ij~ is also localized in the other boundary (s, t) = (N s , N s ). 

The finite terms and linearly divergent terms should be calculated by repeating the 
numerical integration 0(N 2 ) times. However, as can be seen in the next section, such huge 
number of integrations can be avoided for the wave-function renormalization of the quark 
field. On the other-hand, since the structures of I^ nite and M± are useful to understand the 
domain-wall QCD more deeply, we will give them in the separated paper. 



IV. RENORMALIZATION OF QUARKS FIELD 

The result obtained in the previous section is summarized in the following form of the 
effective action for 2 point function with the scale p 2 = (fxa) 2 at 1-loop level: 



r (2) 



where 



H-P)s [il&v (z + p + + z-p_) + W + P + + w P-] st i>(p) 



Z± = 1 + 9^(1^ + 1^ + 1^) 
W ± = W ± (0) + g 2 C 2 {M taA + M±) 



(41) 

(42) 
(43) 



with 



— T5 st = 
2 ' 

-2Ti ( = • 



-0.0773065 



s,t 



-0.3092245 



s,t- 



(44) 
(45) 



^tad(M) = 

M tad (s,t) = 

The expressions for , /^ ite and M ± are given in the previous section. In this section 
we consider the renormalization of zero modes, which is interpolated by the quark field: 
q(p) = P+ip(p)i + P-ip(p)N a ■ Here we only present the results and give the detail of 
derivations in appendix C. 
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A. Diagonalization of mass matrix and stability of zero modes 

For the renormalization of zero modes, it is better to use new basis, ip d {p) which diago- 
nalize mass matrices W . This basis are given by the relation that 

^(p) = U s , t P + 4>(p) + V s , t P^(p), (46) 

where unitary matrices U and V satisfy 

UW~W + U Y 

vw + w~v ] ~ 



s,t 



In our notation the mass eigen- value squared M 2 is arranged in such a way that M^ s =0, 
and we can take U and V real matrices without loss of generality. 
We calculate U and V at 1-loop level: 

U = {l + ff 2 U 1 )U , V = (l + g 2 V 1 )V , (47) 

where tree level matrices U and V are analytically obtained in the large N s limit as follows: 

_( (2/N s y/*sma s (N s + l-t) s ? N. 

[Uokt - { (i - wiy/ 2 wt 1] s = n s (48) 

and [V ] s>t = [Uo\ s ,N s +i-t , where w = 1 — M with M = M + A(u — 1). Here u — 1 for the 
naive perturbation theory, while u — 1 — g 2 C2T/2 for the tadpole improved perturbation 
theory. Hereafter we will call both cases "the tree level" and will not distinguish the two cases 
unless necessary. If we expand the mass eigen-value squared as (M 2 ) s = (Mq) s + g 2 (M 2 ) s , 
the tree level one is related to the phase factor a s such that 2w coso; s = 1 + Wq — (Mq) 8 . 
This phase factor, which also satisfies sinc^A^ = Wosma s (N s + 1), is explicitly given as 
a s = ns/N s in the large N s limit. It is also shown that U and V diagonalize W ± itself such 
that [V W + UlU = [UoW-VjLt = (M ) s S s , t ._ 

We now consider W , which is denoted as W = W^+g 2 W^, where (W^) Sjt = W ± (0) = 
St,s±i - w and g 2 (Wt) s ,t = g 2 C 2 (M ± + M tad ) s , t + 4(1 - u)5 s , t . To diagonalize W T ■ W ± at 
1-loop order, U\ and V\ should satisfy 

(C/i) S)t (M 2 ), + (M 2 ) s (Ul)s,t 

+ (U Wr Vj ■ V W£ul) Sit + (U W ~Vj ■ V W?Ut) 8 , t = (M 2 ) s 5 s , t , (49) 

(Vi)., t (M?)t + M).(tf).,* 

+ (Wi + C/ f • U WoV^) 8 , t + (VoW+U^, ■ U WrvJ) S:t = (M 2 ) s 5 s , t . (50) 

Using the fact that (Ui,Vi) S;t = —(Ui,Vi) tjS implied by the unitarity and the reality, and 
VqWqUq and UqWqVq are diagonal, we can easily solve the above equation as 

{Ul)s t = W)t(Wi)*, + (M )s(W 1 ) s , 



(Vi) 



s,t 



(M 2 ) s - (Mo 2 ), 
(W^sAMoh + (W 1 ) ttS (M ) t 
(M$) s - (M 2 ) t 
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for s 7^ t, and 

(M 2 ) s = 2{W 1 ) S , S {M Q ) S , {U x ) atS = (Vi),,, = 0, (51) 

where W 1 = WiE/J. Th e mass eigen- value squared M 2 = (M^% + g 2 (M 2 ) s obtained above 
leads to the mass eigen-value M s itself: M s = (M ) s + <7 2 (Wi) SjS . Note that M^ s = 
since (Mo)jv s = and (Wi)n s ,n s = in the large N s limit as is shown in appendix C. This 
result explicitly demonstrates the stability of the zero modes against 1-loop corrections in 
domain-wall QCD. As in the case at the tree level, it is shown that 

(VW + U^) s , t = (UW~V% t = M s 5 s , t + 0{g A ). (52) 



B. Wave function renormalization for quark fields 



After diagonalization of the mass matrix, the effective action for the zero mode field 
^ d {p)N a = Xo(p) becomes 



Xo(-P) [w^ { Z + P + + Z - P -)] Xo(p) 



where 



7- i 2 r T + ? c ' 
Z± = 1 - g C 2 — + 



, |67r |(log7r 2 + ^-logM 2 )+/(/l)v...v. 



(53) 



(54) 



with I d = C 2 (U Ifi nite Uti) and I d = C 2 (V / finitc Vq) . Since the interpolating quark 
field q(p) is expressed as q(p) = (U Nsl P + + V NatNt P-)xo(p), and (xo(p)Xo(-p)) = 



Z + Z_ 



, we obtain 



(q(p)q(-p)) 



U 2 

U N a l 



p_ 



p2 



(55) 



— — t^y J) 

Therefore, the renormalized quark field Q(p), which satisfies (Q(p)Q(—p)) = tt- -, is 

/>- 

given by Q(p) = [(Z^ 2 P + + (Z F ) 1 / 2 P_]g(p) with Z^ = J±- and Z„ = Since 



an explicit evaluation shows that (I. 



N s lj 



2 <^ ^ t 

U N S 1 ^N S ,N S 

+)n s ,n s = (H)n.,n. = I d , thus Z + = Z_ = Z, and 

Z 



N S ,N S 



1 — Wq , we finally obtain Zp = Z F = Zp = ^ where 



Z = l~ g 2 C 2 T - + Yf^OogTr 2 + \- log(/ia) 2 ) + g 2 I d . 



(56) 



Here one unknown constant I d is given by 
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c ( dH 

x - sin 2 Z M 

1 



(iE[ sin2 ^(^ + G L ) + 2cos/, 



> - W{l))G 



R 



(wo - W(l))G R + (£ cos 2 Z„/2 - 2 cos 2 / M /2)G L + ^(sin 2 l v /2)G R 



(I 



2\2 



7l Z ~ V 



(57) 



where 



G L = A 
Gr = A 



G 
G 



e a -W 



1 



with 



A i - < 

2VFsmha 
q sinh a — sinh a 

2u> sinh a (cosh a — cosh a) 

and e _a ° = Wo- The numerical value of I d is given in Table | at several values of M, together 
with the total 1-loop renormalization factor Z\ ( Z = 1 + g 2 Zi ) at fia = 1 and the ratio 

C 2 



of the non-tadpole contribution (-Zi)non-tad = I d 



16n 2 



; (log7T 2 + 0.5) = I d + 0.02355 to the 

total one. Note also that the tadpole contributionjs;ives Z tad = —C 2 T/2 = —0.1031. From 
this table, we see that I d is small and depends on M very weakly: The value I d = -0.01945 
at M = 0.05 monotonically increases (; decreases in the absolute value ) to I d = —0.01222 
at M — 0.95. Furthermore the non-tadpole contribution Z non _ ta( j is relatively small: 4% at 
M = 0.05 and 12% at M — 0.95, so that the tadpole contribution becomes dominant at 
all M. This justifies the use of the tree- level result with the tad-pole improvement. Since 
Z\ ~ 0.1, the one- loop correction to the Z factor is about 10 % at g 2 ~ 1.0. 



V. MEAN FIELD ANALYSIS AT FINITE N s 

As seen in the previous sections, due to the presence of off-diagonal terms in the extra 
dimension, analysis of the 1-loop correction to domain-wall quarks becomes too complicated 
to be easily applied to results of the numerical simulations, which should be performed on 
finite N s . In this section we adopt an approximated but simpler method to analyze the 
effect of 1-loop corrections. We call the method the mean field (MF) analysis since the 
link variable U njfl in the fermion action is simply replaced by the mean field u which is 
independent on n and fi. After this replacement the fermion propagator can be explicitly 
calculated and result is identical to the one given in Appendix B with the replacement such 
that x — > ux and cosp^ — > ucosp^. In perturbation theory this is equivalent to the tree 
level analysis with the tad-pole improvement, which has been shown in the previous section 
to give about 90 % of the wave function renormalization factor at 1-loop level. 
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Since we are interested in the zero mode at s = 1, we set s = t = 1 in the propagator. 

2a of Gl, which is given at non-zero m q by 



In this case the zero mode appears in B e 

(l-We~ a )(l-m 2 ) 



-2a 



2Wsmh{a)F 



(58) 



where 



F = We a - 1 + m\{\ - We~ a ) -Am q -W ■ sinh(a) e - a ^ 
+ e~ 2aNs (l - We~ a + m 2 q (We a - 1)). 

In the small momentum limit, this leads to 

z- 1 



lim B e 



-2a 



X 



(1 — wl) 2 + p 2 uwl 



p 2 + m F (1 



2\2 



+ P 2 U(1 + Wq) 



(59) 



(60) 



where Z~ 



B 



m, 



A 



Au 1 r Au 
1 



and w = 1 - M + 4(1 - it) = 1 - M with 



1-w 2 



1 + m\wl 



iu (l - wl)m 2 + m q WQ s 
x {2iV s (l - wl) -l-w 2 + 2w {l - wl) -N S {1- w 2 ) 2 /w } 



+ wl N °{w' 2 + m 2 - 2N.(1 - wl) - w (l - w 2 ) + N,(l - w 2 ) 2 /w } 



B 



1 



2m q w Ns + w'l 



2N S 



Since the pole in the second factor is in general larger than the physical pole in the first 
factor, we neglect the second factor in the latter analysis. 

Now we use the above formula to understand the behavior of the zero mode observed in 
ref. [[|. For the value of u there are several choices. The tadpole diagram alone gives 



u 



1 - g 2 C 2 T/2 = 1 



0.10307/ ~ exp[-0.10307# 2 ] 



where we may take the bare coupling 2N C / (3 or the renormalized coupling gj-^(Tc/a) for g 2 
in the above formula. Alternatively we may also use the "observed" value of it: it = P 1 / 4 
where P is the average value of the plaquette normalized to unity. We adopt the latter 
one in our analysis. The configurations in ref. || generated at /3 = 5.7 and m q a = 0.01 by 
the dynamical Kogut-Susskind quark action give P = 0.5772, which leads to it = 0.872. In 
ref. two remarkable features are found for the zero mode: no zero mode is observed for 
N s = 4 and the zero mode is observed at M = 1.7 but not at M < 1.0 for N s = 10. To 
explain these we calculate mpasa function of M for both N s = 4 and 10 at m q = 0, 0.01, 
0.02 ,0.03, and plot the results in Fig.^, where solid lines are for N s = 4 and dashed lines 
for N s = 10. Four lines for each iV s correspond to m q =0, 0.01, 0.02, 0.03 from below to 
above around M = 1.5. The result tells us the followings. The allowed range for the light 
fermion is very narrow for N s = 4 (roughly 1.4 < M < 1.6 ). This may be a reason why 
the light state could not be found in the simulation || . Note that the allowed range for the 
zero mode is 0.512 < M < 2.512 in the N s —>■ oo limit. Although the allowed range becomes 
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larger for N s = 10 (1.1 < M < 1.9), no light state appears at M < 1.0, as observed in the 
simulation. Furthermore the order of the fermion mass mp is reversed to the order of the 
current quark mass m q at M < 1.0: m F is largest at m q = 0. The plot also supports the 
fact that the zero mode is observed at M — 1.7 in the simulation. 

As seen in the above the behavior of the numerical simulation is understandable by 
the MF analysis, which can supply useful informations on the tuning of parameters in 
numerical simulations such as N s , M or m q before-hand. For example we may take N s = 4 
for the simulations, which reduces the cost of both CPU time and memory a lot, if M is 
appropriately chosen ( M ~ 1.5 for U = 0.872 ). 

VI. CONCLUSION AND DISCUSSION 

In this paper we calculate one-loop correction to the fermion propagator in the massless 
lattice QCD formulated via domain-wall fermions. We show that the zero mode is stable 
against the one-loop correction: no additative counter term to the quark mass is generated in 
the large N s limit. This property is very different from and superior to the ordinary Wilson 
fermion formulation. We explicitly calculate the wave-function renormalization factor for 
the massless quarks and show that the tadpole contribution becomes dominant at all M. 
We also adopt the mean-field analysis to this model, demonstrating that it can qualitatively 
explain data obtained in the numerical simulation ||. 

Although our results strongly indicate that the domain-wall QCD can avoid the fine 
tuning problem of the quark mass, the mechanism which gives the zero mode in this formu- 
lation has not been fully understood yet. Since our proof for the stability of the zero mode 
contains an explicit calculation at 1-loop (: (Wi)n s ,n s = 0), it can not be easily carried 
over to higher orders. The result of numerical simulation suggests that the zero mode is 
also stable against the non-perturbative dynamics. There may be a yet unknown symmetry 
which ensures the existence of zero mode in the large N s limit. To find such a symmetry 
will be important for our understanding of the formulation 

In this paper only the wave-function renormalization factor is explicitly evaluated. Based 
on the method developed in this paper, it is possible to calculate more complicated quantities 
such as renormalization factors for the quark mass, currents and 4-fermi operators, which 
are necessary to get the continuum physics from numerical simulations. The results of this 
paper also suggest that the smeared quark operator g smear = ^ s (ty*F + ^ s + w s ~ s P^ip s ) may 
give better signals than q = P+ipi + P-tpN a does, since it has a larger overlap to zero modes. 

After this work has been completed, there appears a new paper H] , in which the stability 
of the zero mode is generally considered. 
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APPENDIX A. ACTION AND FEYNMAN RULES 



The gauge part of the action is exactly same as that of the ordinary lattice QCD action 



•Sgauge — E E R e tr (Un,[iUn+0 :fl U n +fi,vU n ^ 



Sfp 



n,fi 



Cn+A^a 1 ( 9M n + \P)) ~ E ab + ^P) ) 4 



5, 



measure 



1 /l-cos(^(n+|/i)ad(r c 

-^LL trln — 7 — " ; — 

1 n n \ [gA c (n + lji)sLd(T c 



(61) 

(62) 
(63) 

(64) 



^ ~z< ' \ > j / ab 

where g is the coupling of the SU(N C ) gauge, (3 = 2N c /g 2 . a is the gauge parameter. The 
actions Spp and ^measure is not needed in our calculation at one loop level. 
The momentum representation of gauge part is 



S 



gauge 



+ S 



GF 



1 

~ 2 
+ ■ ■ 



d A p 



p\ 



1. 

,1 jPfj.Pu 

a 



AM 



(65) 

here + • • • denotes the gluon self interactions which do not come into play in our calculation. 
The fermion-gauge interaction terms in the momentum representation is 

d 4 k d A p d% d% 



S 



int. 



E 

71=1 



(2tt) 4 (2tt) 4 (2tt) 4 (2tt) 



(2vr) 4 5 4 (fc + p + / 1 + ••• + /„ 



x-sTA^ih) ■ ■ ■ A a "(l n )ij(k) s T a > ■ --T a - 



x 



2 V 



_ (Pi*-kp) _|_ 



^{p)s. (66) 



The domain- wall fermion propagator is already given by the Eq. fll2|) . 

The fermion gluon interaction vertices are given by flB5|). Although there are infinite 
number of interactions in the lattice perturbation theory, only two of them are needed for 
the present purpose. One of them is the fermion interaction vertex with one gluon field, 
which is given by 



1 1 

Vi(k,p;l,a;ii) = -igT a {j^ cos -(-k/j, +p fl ) - ir sin -(-k tl + p^)}. 

The other is the vertex with two gluon fields, given by 



(67) 



V 2 (k,p;h,a,l 2 ,b;fJ.) = -g 2 -{T a , T b }{r>V sin -(-k^ + - rcos -(-k^ + p^)}S fiU . (68) 



2" 2 

The gluon propagator is given by 

1 

P 2 



We set a = 1 in this paper. 



Sau ~ (1 - Oi) 



PjjPv 
p 2 



Jab- 



(69) 
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APPENDIX B. DERIVATION OF FREE FERMION PROPAGATOR 



In this appendix we derive the free fermion propagator, used in the text. For the later 
use in perturbative analyses of this model, non-zero current quark mass m q for finite N s is 
considered. See also Refs. . We also derive the propagator with Majorana mass terms, 

which becomes important for the lattice definition of the N — 1 supersymmetric model via 
domain- wall fermions [1C,11]. 

A. Propagator with non-zero m q 

The free fermion propagator has the following form: 

S F (p)s,t = [HTaJV + W-)G R {s,t)P+ + + W+)G L (s,t)P- 



St 



where 



G R {s,t) = ( 1 ) and G L (s,t) ^ 



\f 1 "m m/st \" m ' m / st 

with 

(W£)-,t = (W + ) s , t + m q 5 s , Ns 8 ttl and (W-)., t = + m q 8 sA 8 tiNs (70) 

We first consider G R . The following equation is satisfied for G R . 

£ [(x + W + W~) Sjt + m q (W+6 m . + $s,N s Wr t ) + m 2 q 8 SjNs 8 tNs ] G R (t, u) = 5 SU (71) 
t 

with x = p 2 . Therefore, except s = N s or 1, this equation is satisfied by 

G R (s, t) = G{s, t) + A ++ e a{s+t) + A + _e a ^ + A_ + e a{ ~ s+t) + A__e a( - s -' } (72) 

where 

G(s, t)=A ( e a ^-l s -*l) + e -a(*.-l«-*D) (73) 

becomes a special solution to the equation (x + W + W~)G R = 1, with 

1 + W 2 + x 



cosh a = 



2W{p) 



A = - - (74) 

2Wsinha2sinh(aiV s ) , v ; 

and other terms are general solutions to the equation (x + W + W~)G R = 0. We can fix their 
coefficients A±± by a boundary condition at s = 1: 

(x + W 2 + 1)G R (1, t)-W- G R (2, t)-W-m q - G R (N S , t) = 6 U , (75) 

which is simplified to 
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G R (0,t)-m q G R (N s ,t)=0, (76) 
and another boundary condition at s = N s : 

(x + W 2 )G R (N S , t)-W- G R (N S -l, t )-W-m q - G R (1, t) + m 2 q G R (N s , t) = S N „ t , (77) 
which is reduced to 

G R (N S , t)-W- G R (N S + l,t) + W- m q G R (l, t) - m 2 q G R (N s , t) = 0. (78) 
Plugging eq.(|72|) into eqs.([7(|) and (|78|) leads to 

/ 1 - m q e aNs 1 - m q e- aNs \ ( A ++ A + _ \ 

\ e aNs (l - We a -m 2 q + Wm q e a(1 - N ^) e - aNs (l - We~ a - m\ + Wm q e a{ - N ^) ) \ A_ + A_- J 

_ A (e-" N °-m q e aN °-m q \ 

\1- We~ a -m 2 q + Wm q e- a( - Ns+ V (1 - We a - m 2 + Wm q e a( - N ° + V) J' K ' 

Solving this we obtain 



A ++ \ ( e ~ 2aNs " - We~ a )(l - m 2 q ) 

A_ + ) ~ F\2W sinh(a) (1 - 2m q cosh(aiV s ) + m; 

A+_ \ A f 2W sinh(a)(l - 2m q cosh(aN s ) + m; 



q' 



\A~.) F \(l-e 2aN °)(l-We a )(l-m 2 q ) 
where 

F = e aNa [1 - We a + m 2 q (We- a - 1)] + AWm q sinh(a) + e~ aN ° [We~ a - 1 + m 2 q (l - We a )]. 

(80) 

Similarly, plugging the general solution for Gl 

G L {s, t) = G{s, t) + B ++ e a(s+t) + B + ^e a{s - t] + B_ + e a ^ s+t ^ + B^e^ 3 ^ (81) 
into the boundary conditions 

G L (N B + l,t)-m q G L (l,t) = (82) 
G L (1, t)-W -G L (0,t) + W ■ m q G L (N s , t) - m 2 q G L (l, t) = 0, (83) 

we finally obtain 

A ( ( e - 2aN ° - l)e- a (e- a - W)(l - w 2 ) \ 







\B- +) 













F \ 2iysinh(a)(l - 2m q cosh.(aN s ) + m 
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B. Propagator with the Majorana mass term at N s 



For an application of the free fermion propagator obtained in the domain-wall model, we 
consider the model with the Majorana mass term on the anti-boundary at s — N s , which 
has been proposed for a lattice definition of the N — 1 super Yang- Mills theory (T^JD! • Here 



we set m q = 0. 



A free fermion action of the model with the Majorana mass uiq can be written in the 
momentum space as 



S = h(-p) s D Sjt {p)*(p) 



where 



and 



D{p) = T (p) + m X 



Do(p) 






-A)(-p) 



T + m 5 2 P + IP. 



with (5 2 ) Stt = 5 s , Ns 8 Ns<t , and 



P+ 
P_ 



P_ 



P_ 
P + 



I 2 

h o 



( o 2 \ 

o 2 

a 2 

V a 2 J 



in terms of 8 x 8 matrices. Here 

A>(p) = *7mP„ + W + P + + W-P-, 

is an inverse of the massless free fermion propagator in the domain- wall QCD. 
By expanding D~ l in mo and rearranging it we obtain 

oo 

IT 1 = Yti-T^moP+IP^rT, 1 

n=0 
oo 



n=0 



where Z = IP_5T Q l 5P+. 

Using Z 2 = — x(G n(p) N at N s ) 2 P + and suming over n, we finally get 



D- 1 = T, 1 



-moTo^P+IP-dTu 1 + mlT^dP+ZIP.STo 1 



x 



(84) 



(85) 



(86) 



(87) 



(88) 



1 + m§x(Gfl(p) A r S)JVs 
Explicitly this formula gives, in terms of 2 x 2 block notations, 



|2 ' 



39) 
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D(p)ii = ~D{-p) 22 l = &(pM-p)) 

= (-i 7l ji)[Z + (p)P + + Z_(p)P_] + M + (p)P + + M_P_ 



(90) 



where 



Z+(p) s ,t = 




Z -(P)-S,t = 


G L (p)s,t + J 


M+{p) s ,t = 


{W-G R (p) a , 


M-(p) 8tt = 


(W + G L (p) s , t 



I ° 2 7^ G R (p) StNs G R (p) 



N s ,t 



UlnG 



-(W-G R (p)) s>Ns (G R (p)W + ) Ns 



™ Y r2 (W-G R (p)) StNa G R (p) Nsit 
1 + ijiqxG 2 

m^xG 

1 + rrinxG 2 



G r {p) s ,nXGr{p)W" 



IN. A 



with G = G r (p)n s ,n s - Similarly 



(£»(p)- ] 
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(tp(p)si>(-p)t) 

m 



1 + mniG 2 



xG R (p) StNs G R (p) Ns , t I 2 P_ - ip„l^G R (p) SiNs (G R (p)W- 



)N a ,thP+ 



+ipM W ~ G M)s,N s G R (p) Ns , t I 2 P^ + (^-G fl (p)) s ^ s (G fl (p)^ + )^/ 2 P + 



(91) 



and 

(D( P )21 1 )s,t 



(1/)(p),1/)(-p)t) 

1 + mlxG 2 



xG R (p) S!Na G R (p) Ns!t I 2 P+ + ipX G R^)^AG R (p)W + ) Ns ^I 2 P. 



fy%(W-G R (p)) atN ,G R (p) N . )t I 2 P + + (W-G R (p)) s , Na (G R (p)W + ) Ns , t I 2 P. 



• (92) 



See ref. 11 for an application of this result. 



APPENDIX C. PROPERTIES OF DIAGONALIZATION MATRICES 

In this appendix we derive several properties of diagonalization matrices, U and V, which 
are used for the renormalization of quarks fields. 

Let us consider the tree level diagonalization of matrices (W^ ■ W^ 1 ). To diagonalize 
(W^ ■ Wq) } we have to solve the eigen-value problems (W T • W^) 8 j<l)±(t) = (M^^s), 
then Uq and V are given by normalized eigenvectors 0±: (Uo) s>t = 4> s + {t) and (V ) s>t = (fi s _{t). 
The two eigen-state equations lead to the same equation 



- w + 1) + ^(s - 1)) + (l + wl - (M 2 ) t ) = 

but with different boundary conditions: 

- w^UO) + <p + (l) = 0, <pUN. + 1) = 



(93) 



(94) 
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or 



-w <i> i _(N a + l) + <j>4N s )=0, l _(O)=O. (95) 

Therefore, once <fr l + (s) is known, the other is easily obtained through 0!_(s) = 0^(iV s + 1 — s). 
Hereafter we consider (f> l + (s) only and drop the suffices + and i. 

There are two types of solutions to the eigen-state equation. For (Mg)j < (1 — w ) 2 we 

1 _|_ w 2 _ (M 2 )i 

have a damping solution (j>(s) = Ae~ as with cosh a = — . The first boundary 

condition leads to e~ a = wq. This implies (Mq) = 0, and therefore 4>(s) is nothing but the 
zero mode solution of the domain-wall QCD. For this solution w should satisfy Wq < 1 ( 
< M < 2). The other boundary condition can be satisfied in the large N s limit. The 
normalization constant becomes A = (1 — Wq) 1 ^ 2 . Note that there are no other damping 
solutions which satisfy the first boundary condition. 

If the eigen-value is in the region (1 — wq) 2 < (M 2 )j < (1 + w ) 2 , we have an oscillating 

1 _|_ w 2 _ (M 2 )i 

solution 4>(s) = Ae ias + Be~ tas with cos a = — . The two boundary conditions 

2w 

imply 



e ta — Wq e la — Wq 
e ia(N a +l) e -ia(N s +l) 




0. (96) 



The existence of the non-trivial solution requires Wosina(N s + 1) = sinaiV s , which leads 
to 0(s) = -Ae ia{Ns+1) sina(A?" s + 1 - s) = A sina(iV s + 1 - s). Without loss of generality 
we can take real A , and the normalization condition gives A = (2/N s ) 1 ^ 2 (l + 0(1 /N s )). 
Setting a = a/N s we reduce the equation for a to w sina = sin a in the large N s limit. 
The solutions a = nn with integer n to this equation is translated to iV s — 1 independent 
solutions: a = rcn/N s with n — 1, 2, • • • N s — 1. (Note that < a s < n since sina s > 0.) 
Therefore, all eigen-values and eigen-vectors are now obtained, giving 

/ (2/N s y/ 2 sma s (N s + l-t) s ? N s 

[f/o ^ = \ (i-w 2 y/ 2 wt 1] s = n s > (97) 

and [Vo] a>t = [U ] s ,N s +i-t- 

Next we prove some properties of U and V . It is noted that U and V can also 
diagonalize Wq : 

(VoW+Vt) =6,, t f, (98) 



where 



, _ f w cosa s (N s + 1) - cosa s A^ s s ^ N s , , 

Is ~ \ s = N s ' iyyj 

Using the equation for a s (s ^ N s ) we can show 

f 2 = vol + 1 — 2wq [sin a s N s sin a(N s + 1) + cosa s A^ s coso; s (A^ s + 1)] 

= wl + l-2w cosa s = (M 2 ) s . (100) 



19 



This proves f s = (M ) s for all s. 

It is also important to note that U (V ) diagonalizes l£ (iw) terms, since 



(U vtf*- 2 ui) at = 5 s>t 5 s , Ns (l - O E^o (s+4 ~ 2) + 0(1/N S ) 

s,t 

= S s , t S s , Ns (l - w 2 )' 1 + 0{l/N s ). (101) 



Now let us consider quantities including g 2 contributions. We first show that U and V 
diagonalize W at this order: 

V ■ W + ■ = (1 + g 2 V 1 )V (W+ + g 2 W+)Ul(l + g 2 U{) 

= v ■ w+ul + g 2 |y x v w+ui + v w u}ul + (102) 

where the coefficient of the g 2 term is simplified to 

(Yi).AMo) t + (M )s(ul) st + {wfu = (wf) s , t (i + { -^0^f) = ( 103 ) 

for s 7^ t, and becomes (W% ) S)S for s = t. Eq. (|102|) then becomes 

= (m + 9 2 Wx + } 1. 

It is necessary for the stability of the zero mode to show that (Wi + )n Si n s = 0. This can be 
proven as follows. 

(m + ) Ns , Ns = (1 - wl) £ «^'- W).,^ 1 (104) 

where Wf 1- is composed of the sum Gl, Gr, WqGl and W ~Gr but 

J^-'GfotK -1 = 0(N.w»') - 

s,t 

for all G = Gl, Gr, WqGl and W ~Gr in the large iV s limit. Therefore eq.( |104| ) vanishes. 
For the wave-function renormalization factor we have to know 

U Ns ,i = V Ns>Na = {U ) Nsjl + g 2 E {Ui)N.,t(Uo)t,i- 

t^N s 

Fortunately, since (Z7o)t,i = (2/iV,,) 1 / 2 sin a t (l — 1) = 0, there is no order g 2 contribution and 
it becomes Un b ,i = (1 - Wq) 1/2 - 

Finally we would like to evaluate I± from I finite . If we define 

(F(s,t))u = J2( U o)N s , s F(s,t)(U ) Ns , 

s,t 

and 

(F( S ,t)) v = J2(V ) Ns ,sF(s,t)%)N s ,t, 

s.t 
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we can show 



(e- als -%,v = (1 



W r . 



sinh ao — sinh a 



2wq sinh ct (cosh a — cosh a) 



with e" 



w , 



-a(s+t-2) 



-a(2JV s -s-t) 



)y = (l-^ 2 )- 



,2a 



e a - w ) 2 



and 



Using these formula we obtain 

sinh ao — sinh a e a — W 



(G L (s,t))u = 



2W sinh a 



2w sinh a (cosh a — cosh a) e 01 — W (e a — w ) 2 



= (G R (s,t)) v = G L 



(G R (s,t))u 



l-w 2 
2W sinh a 



sinh a — sinh a 



2w sinh a (cosh a — cosh a) (e a — wq)' 



= (G L (s,t)) v = G R 



(105) 



(106) 



and 



(W+G L (s,t)) u>v = (W -G R (s,t)) uy = (w - W)G R . 

The explicit expression for I^ nite is reduced to the final result in terms of G L / R : 1^ = 1^ = I d 
where 



I d = C 2 



dH f 1 



X) [sin 2 l„(G R + G£) + 2 cos Z M (w - W(Z))G fl 



+ sin 2 / M 



V 2(/ 2 ) 2 
1 



(2tt) 4 1 8/2 ^ 

(w - W(1))Gr + (]T cos 2 Z„/2 - 2 cos 2 IJ2)(T L + ]T(sin 2 Z„/2)Gr 



6(n 2 -l 2 ) 



(107) 
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FIGURES 




j^half — circle 

<- 1 



FIG. 1. Diagrams which contribute to the one-loop correction to the fermion propagator. 
Above: Tadpole diagram. Below: Half-circle diagram. 
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FIG. 2. The fermion mass mp obtained in the mean- field approximation as a function of M 
for N s =4 (solid lines) and N S =W (dashed line), at m q = 0, 0.01, 0.02, 0.03 from below to above 
around M = 1.5. 
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TABLES 



TABLE I. Value of I d vs M, together with Z x and {Z 1 ) non _ ud /Z 1 



M 




Zi 


l^ljnon-tad/^l 


0.05 


-0.01945(5) 


-0.09897 


0.041 


0.10 


-0.01871(5) 


-0.09822 


0.049 


0.15 


-0.01804(5) 


-0.09756 


0.056 


0.20 


-0.01744(5) 


-0.09696 


0.063 


0.25 


-0.01688(5) 


-0.09640 


0.069 


0.30 


-0.01636(5) 


-0.09589 


0.075 


0.35 


-0.01588(5) 


-0.09541 


0.080 


0.40 


-0.01544(5) 


-0.09496 


0.085 


0.45 


-0.01502(5) 


-0.09454 


0.090 


0.50 


-0.01463(5) 


-0.09415 


0.095 


0.55 


-0.01426(5) 


-0.09378 


0.099 


0.60 


-0.01392(5) 


-0.09345 


0.103 


0.65 


-0.01361(5) 


-0.09313 


0.107 


0.70 


-0.01332(5) 


-0.09284 


0.110 


0.75 


-0.01305(5) 


-0.09257 


0.113 


0.80 


-0.01281(5) 


-0.09233 


0.116 


0.85 


-0.01259(5) 


-0.09211 


0.119 


0.90 


-0.01239(5) 


-0.09191 


0.121 


0.95 


-0.01222(5) 


-0.09174 


0.124 
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